Persistent DNA Methylation Changes across the First Year of Life and Prenatal NO2 Exposure in a Canadian Prospective Birth Study.

BACKGROUND: Evidence suggests that prenatal air pollution exposure alters DNA methylation (DNAm), which could go on to affect long-term health. It remains unclear whether DNAm alterations present at birth persist through early life. Identifying persistent DNAm changes would provide greater insight into the molecular mechanisms contributing to the association of prenatal air pollution exposure with atopic diseases. OBJECTIVES: This study investigated DNAm differences associated with prenatal nitrogen dioxide (NO2) exposure (a surrogate measure of traffic-related air pollution) at birth and 1 y of age and examined their role in atopic disease. We focused on regions showing persistent DNAm differences from birth to 1 y of age and regions uniquely associated with postnatal NO2 exposure. METHODS: Microarrays measured DNAm at birth and at 1 y of age for an atopy-enriched subset of Canadian Health Infant Longitudinal Development (CHILD) study participants. Individual and regional DNAm differences associated with prenatal NO2 (n=128) were identified, and their persistence at age 1 y were investigated using linear mixed effects models (n=124). Postnatal-specific DNAm differences (n=125) were isolated, and their association with NO2 in the first year of life was examined. Causal mediation investigated whether DNAm differences mediated associations between NO2 and age 1 y atopy or wheeze. Analyses were repeated using biological sex-stratified data. RESULTS: At birth (n=128), 18 regions of DNAm were associated with NO2, with several annotated to HOX genes. Some of these regions were specifically identified in males (n=73), but not females (n=55). The effect of prenatal NO2 across CpGs within altered regions persisted at 1 y of age. No significant mediation effects were identified. Sex-stratified analyses identified postnatal-specific DNAm alterations. DISCUSSION: Regional cord blood DNAm differences associated with prenatal NO2 persisted through at least the first year of life in CHILD participants. Some differences may represent sex-specific alterations, but replication in larger cohorts is needed. The early postnatal period remained a sensitive window to DNAm perturbations. https://doi.org/10.1289/EHP13034.

Epigenetic clocks, sex markers and age-class diagnostics in three harvested large mammals.

The development of epigenetic clocks, or the DNA methylation-based inference of age, is an emerging tool for ageing in free ranging populations. In this study, we developed epigenetic clocks for three species of large mammals that are the focus of extensive management throughout their range in North America: white-tailed deer, black bear and mountain goat. We quantified differential DNA methylation patterns at over 30,000 cytosine-guanine sites (CpGs) from tissue samples of all three species (black bear n = 49; white-tailed deer n = 47; mountain goat n = 45). We used a penalized regression model (elastic net) to build explanatory (black bear r = .95; white-tailed deer r = .99; mountain goat r = .97) and robust (black bear Median Absolute Error or MAE = 1.33; white-tailed deer MAE = 0.29; mountain goat MAE = 0.61) models of age or clocks. We also characterized individual CpG sites within each species that demonstrated clear differences in methylation levels between age classes and sex, which can be used to develop a suite of accessible diagnostic markers. This tool has the potential to contribute to wildlife monitoring by providing easily obtainable representations of age structure in managed populations.

Novel DNA methylation changes in mouse lungs associated with chronic smoking.

Smoking is a potent cause of asthma exacerbations, chronic obstructive pulmonary disease (COPD) and many other health defects, and changes in DNA methylation (DNAm) have been identified as a potential link between smoking and these health outcomes. However, most studies of smoking and DNAm have been done using blood and other easily accessible tissues in humans, while evidence from more directly affected tissues such as the lungs is lacking. Here, we identified DNAm patterns in the lungs that are altered by smoking. We used an established mouse model to measure the effects of chronic smoke exposure first on lung phenotype immediately after smoking and then after a period of smoking cessation. Next, we determined whether our mouse model recapitulates previous DNAm patterns observed in smoking humans, specifically measuring DNAm at a candidate gene responsive to cigarette smoke, Cyp1a1. Finally, we carried out epigenome-wide DNAm analyses using the newly released Illumina mouse methylation microarrays. Our results recapitulate some of the phenotypes and DNAm patterns observed in human studies but reveal 32 differentially methylated genes specific to the lungs which have not been previously associated with smoking. The affected genes are associated with nicotine dependency, tumorigenesis and metastasis, immune cell dysfunction, lung function decline, and COPD. This research emphasizes the need to study CS-mediated DNAm signatures in directly affected tissues like the lungs, to fully understand mechanisms underlying CS-mediated health outcomes.

DNA methylation is not associated with sensitization to or dietary introduction of highly allergenic foods in a subset of the CHILD cohort at age 1 year.

Background: In the first year of life, DNA methylation (DNAm) patterns are established and are particularly susceptible to exposure-induced changes. Some of these changes may leave lasting effects by persistently altering gene expression or cell type composition or function, contributing to disease. Objectives: In this discovery study, we investigated DNAm associations with sensitization to peanut, egg, or cow's milk and hypothesized that genes demonstrating DNAm differences in immune cells may play a role in the development of food sensitization. Methods: Infant sensitization (a skin prick test wheal size that is at least 2 mm greater than the negative control) was measured to peanut, egg, and cow's milk at age 1 year, and ages of food introduction were reported prospectively. PBMC DNAm was measured in blood samples at 1 year in 144 infants, oversampled for atopy or wheeze. Statistical analysis of Illumina 450k array DNAm data was conducted in R with adjustment for clinical and genetic covariables and a minimum effect size of 1%, false discovery rate of 5%, and medium-confidence false discovery rate threshold of 20%. Results: There were no DNAm differences between infants with and without peanut, egg, or cow's milk sensitization. Borderline significant sites with high effect sizes were enriched for methylation quantitative trait loci, hinting at genetic factors influencing DNAm at these sites. DNAm patterns did not differ by peanut or egg introduction before or after 12 months. Conclusion: This small pilot study did not show differences in methylation by food sensitization or introduction, but it did demonstrate DNAm patterns linked to genetic variants.

Early-life exposure to cigarette smoke primes lung function and DNA methylation changes at Cyp1a1 upon exposure later in life.

Prenatal and early-life exposure to cigarette smoke (CS) has repeatedly been shown to induce stable, long-term changes in DNA methylation (DNAm) in offspring. It has been hypothesized that these changes might be functionally related to the known outcomes of prenatal and early-life CS exposure, which include impaired lung development, altered lung function, and increased risk of asthma and wheeze. However, to date, few studies have examined DNAm changes induced by prenatal CS in tissues of the lung, and even fewer have attempted to examine the specific influences of prenatal versus early postnatal exposures. Here, we have established a mouse model of CS exposure which isolates the effects of prenatal and early postnatal CS exposures in early life. We have used this model to measure the effects of prenatal and/or postnatal CS exposures on lung function and immune cell infiltration as well as DNAm and expression of Cyp1a1, a candidate gene previously observed to demonstrate DNAm differences on CS exposure in humans. Our study revealed that exposure to CS prenatally and in the early postnatal period causes long-lasting differences in offspring lung function, gene expression, and lung Cyp1a1 DNAm, which wane over time but are reestablished on reexposure to CS in adulthood. This study creates a testable mouse model that can be used to investigate the effects of prenatal and early postnatal CS exposures and will contribute to the design of intervention strategies to mediate these detrimental effects.NEW & NOTEWORTHY Here, we isolated effects of prenatal from early postnatal cigarette smoke and showed that exposure to cigarette smoke early in life causes changes in offspring DNA methylation at Cyp1a1 that last through early adulthood but not into late adulthood. We also showed that smoking in adulthood reestablished these DNA methylation patterns at Cyp1a1, suggesting that a mechanism other than DNA methylation results in long-term memory associated with early-life cigarette smoke exposures at this gene.

Integrative analysis reveals novel associations between DNA methylation and the serum metabolome of adolescents with type 2 diabetes: A cross-sectional study.

Objective: Rates of type 2 diabetes (T2D) among adolescents are on the rise. Epigenetic changes could be associated with the metabolic alterations in adolescents with T2D. Methods: We performed a cross sectional integrated analysis of DNA methylation data from peripheral blood mononuclear cells with serum metabolomic data from First Nation adolescents with T2D and controls participating in the Improving Renal Complications in Adolescents with type 2 diabetes through Research (iCARE) cohort study, to explore the molecular changes in adolescents with T2D. Results: Our analysis showed that 43 serum metabolites and 36 differentially methylated regions (DMR) were associated with T2D. Several DMRs were located near the transcriptional start site of genes with established roles in metabolic disease and associated with altered serum metabolites (e.g. glucose, leucine, and gamma-glutamylisoleucine). These included the free fatty acid receptor-1 (FFAR1), upstream transcription factor-2 (USF2), and tumor necrosis factor-related protein-9 (C1QTNF9), among others. Conclusions: We identified DMRs and metabolites that merit further investigation to determine their significance in controlling gene expression and metabolism which could define T2D risk in adolescents.

In-utero exposure to indoor air pollution or tobacco smoke and cognitive development in a South African birth cohort study.

BACKGROUND AND AIMS: There is increasing evidence indicating that air pollution exposure is associated with neuronal damage. Since pregnancy is a critical window of vulnerability, air pollution exposure during this period could have adverse effects on neurodevelopment. This study aims 1) to analyze associations of prenatal exposure to indoor air pollution (particulate matter with diameters ≤10 µm, PM10) and tobacco smoke with neurodevelopment and 2) to determine whether these associations are mediated by deviations of epigenetic gestational age from chronological gestational age (ΔGA). METHODS: Data of 734 children from the South African Drakenstein Child Health Study were analyzed. Prenatal PM10 exposure was measured using devices placed in the families' homes. Maternal smoking during pregnancy was determined by maternal urine cotinine measures. The Bayley Scales of Infant and Toddler Development III (BSID-III) was used to measure cognition, language and motor development and adaptive behavior at two years of age. Linear regression models adjusted for maternal age, gestational age, sex of child, ancestry, birth weight/length, and socioeconomic status were used to explore associations between air pollutants and BSID-III scores. A mediation analysis was conducted to analyze if these associations were mediated by ΔGA using DNA methylation measurements from cord blood. RESULTS: An increase of one interquartile range in natural-log transformed PM10 (lnPM10; 1.58 µg/m3) was significantly associated with lower composite scores in cognition, language, and adaptive behavior sub-scores (composite score ß-estimate [95%-confidence interval]: -0.950 [-1.821, -0.120]). Maternal smoking was significantly associated with lower adaptive behavior scores (-3.386 [-5.632, -1.139]). Associations were not significantly mediated by ΔGA (e.g., for PM10 and cognition, proportion mediated [p-value]: 4% [0.52]). CONCLUSION: We found an association of prenatal exposure to indoor air pollution (PM10) and tobacco smoke on neurodevelopment at two years of age, particularly cognition, language, and adaptive behavior. Further research is needed to understand underlying biological mediators.

CEBPß regulation of endogenous IGF-1 in adult sensory neurons can be mobilized to overcome diabetes-induced deficits in bioenergetics and axonal outgrowth.

Aberrant insulin-like growth factor 1 (IGF-1) signaling has been proposed as a contributing factor to the development of neurodegenerative disorders including diabetic neuropathy, and delivery of exogenous IGF-1 has been explored as a treatment for Alzheimer's disease and amyotrophic lateral sclerosis. However, the role of autocrine/paracrine IGF-1 in neuroprotection has not been well established. We therefore used in vitro cell culture systems and animal models of diabetic neuropathy to characterize endogenous IGF-1 in sensory neurons and determine the factors regulating IGF-1 expression and/or affecting neuronal health. Single-cell RNA sequencing (scRNA-Seq) and in situ hybridization analyses revealed high expression of endogenous IGF-1 in non-peptidergic neurons and satellite glial cells (SGCs) of dorsal root ganglia (DRG). Brain cortex and DRG had higher IGF-1 gene expression than sciatic nerve. Bidirectional transport of IGF-1 along sensory nerves was observed. Despite no difference in IGF-1 receptor levels, IGF-1 gene expression was significantly (P < 0.05) reduced in liver and DRG from streptozotocin (STZ)-induced type 1 diabetic rats, Zucker diabetic fatty (ZDF) rats, mice on a high-fat/ high-sugar diet and db/db type 2 diabetic mice. Hyperglycemia suppressed IGF-1 gene expression in cultured DRG neurons and this was reversed by exogenous IGF-1 or the aldose reductase inhibitor sorbinil. Transcription factors, such as NFAT1 and CEBPß, were also less enriched at the IGF-1 promoter in DRG from diabetic rats vs control rats. CEBPß overexpression promoted neurite outgrowth and mitochondrial respiration, both of which were blunted by knocking down or blocking IGF-1. Suppression of endogenous IGF-1 in diabetes may contribute to neuropathy and its upregulation at the transcriptional level by CEBPß can be a promising therapeutic approach.

Immune cell type and DNA methylation vary with reproductive status in women: possible pathways for costs of reproduction.

BACKGROUND: Consistent with evolutionarily theorized costs of reproduction (CoR), reproductive history in women is associated with life expectancy and susceptibility to certain cancers, autoimmune disorders and metabolic disease. Immunological changes originating during reproduction may help explain some of these relationships. METHODOLOGY: To explore the potential role of the immune system in female CoR, we characterized leukocyte composition and regulatory processes using DNA methylation (DNAm) in a cross-sectional cohort of young (20-22 years old) women differing in reproductive status. RESULTS: Compared to nulliparity, pregnancy was characterized by differential methylation at 828 sites, 96% of which were hypomethylated and enriched for genes associated with T-cell activation, innate immunity, pre-eclampsia and neoplasia. Breastfeeding was associated with differential methylation at 1107 sites (71% hypermethylated), enriched for genes involved in metabolism, immune self-recognition and neurogenesis. There were no significant differences in DNAm between nulliparous and parous women. However, compared to nullipara, pregnant women had lower proportions of B, CD4T, CD8T and natural killer (NK) cells, and higher proportions of granulocytes and monocytes. Monocyte counts were lower and NK counts higher among breastfeeding women, and remained so among parous women. IMPLICATIONS: Our findings point to widespread differences in DNAm during pregnancy and lactation. These effects appear largely transient, but may accumulate with gravidity become detectable as women age. Nulliparous and parous women differed in leukocyte composition, consistent with more persistent effects of reproduction on cell type. These findings support transient (leukocyte DNAm) and persistent (cell composition) changes associated with reproduction in women, illuminating potential pathways contributing to CoR. Lay Summary: Evolutionary theory and epidemiology support costs of reproduction (CoR) to women's health that may involve changes in immune function. We report differences in immune cell composition and gene regulation during pregnancy and breastfeeding. While many of these differences appear transient, immune cell composition may remain, suggesting mechanisms for female CoR.

Newborn differential DNA methylation and subcortical brain volumes as early signs of severe neurodevelopmental delay in a South African Birth Cohort Study.

OBJECTIVES: Early detection of neurodevelopmental delay is crucial for intervention and treatment strategies. We analysed associations between newborn DNA methylation (DNAm), neonatal magnetic resonance imaging (MRI) neuroimaging data, and neurodevelopment. METHODS: Neurodevelopment was assessed in 161 children from the South African Drakenstein Child Health Study at 2 years of age using the Bayley Scales of Infant and Toddler Development III. We performed an epigenome-wide association study of neurodevelopmental delay using DNAm from cord blood. Subsequently, we analysed if associations between DNAm and neurodevelopmental delay were mediated by altered neonatal brain volumes (subset of 51 children). RESULTS: Differential DNAm at SPTBN4 (cg26971411, Δbeta = -0.024, p-value = 3.28 × 10-08), and two intergenic regions (chromosome 11: cg00490349, Δbeta = -0.036, p-value = 3.02 × 10-08; chromosome 17: cg15660740, Δbeta = -0.078, p-value = 6.49 × 10-08) were significantly associated with severe neurodevelopmental delay. While these associations were not mediated by neonatal brain volume, neonatal caudate volumes were independently associated with neurodevelopmental delay, particularly in language (Δcaudate volume = 165.30 mm3, p = 0.0443) and motor (Δcaudate volume = 365.36 mm3, p-value = 0.0082) domains. CONCLUSIONS: Differential DNAm from cord blood and increased neonatal caudate volumes were independently associated with severe neurodevelopmental delay at 2 years of age. These findings suggest that neurobiological signals for severe developmental delay may be detectable in very early life.

Maternal psychosocial risk factors and child gestational epigenetic age in a South African birth cohort study.

Accelerated epigenetic aging relative to chronological age has been found to be associated with higher risk of mortality in adults. However, little is known about whether and how in utero exposures might shape child gestational epigenetic age (EA) at birth. We aimed to explore associations between maternal psychosocial risk factors and deviation in child gestational EA at birth (i.e., greater or lower EA relative to chronological age) in a South African birth cohort study-the Drakenstein Child Health Study. Maternal psychosocial risk factors included trauma/stressor exposure; posttraumatic stress disorder (PTSD); depression; psychological distress; and alcohol/tobacco use. Child gestational EA at birth was calculated using an epigenetic clock previously devised for neonates; and gestational EA deviation was calculated as the residuals of the linear model between EA and chronological gestational age. Bivariate linear regression was then used to explore unadjusted associations between maternal/child risk factors and child gestational EA residuals at birth. Thereafter, a multivariable regression method was used to determine adjusted associations. Data from 271 maternal-child dyads were included in the current analysis. In the multivariable regression model, maternal PTSD was significantly and negatively associated with child gestational EA residuals at birth (ß = -1.95; p = 0.018), controlling for study site, sex of the child, head circumference at birth, birthweight, mode of delivery, maternal estimated household income, body mass index (BMI) at enrolment, HIV status, anaemia, psychological distress, and prenatal tobacco or alcohol use. Given the novelty of this preliminary finding, and its potential translational relevance, further studies to delineate underlying biological pathways and to explore clinical implications of EA deviation are warranted.

Methylation of cysteinyl leukotriene receptor 1 genes associates with lung function in asthmatics exposed to traffic-related air pollution.

Air pollution is associated with early declines in lung function and increased levels of asthma-related cysteinyl leukotrienes (CysLT) but a biological pathway linking this rapid response has not been delineated. In this randomized controlled diesel exhaust (DE) challenge study of 16 adult asthmatics, increased exposure-attributable urinary leukotriene E4 (uLTE4, a biomarker of cysteinyl leukotriene production) was correlated (p = 0.04) with declines in forced expiratory volume in 1-second (FEV1) within 6 hours of exposure. Exposure-attributable uLTE4 increases were correlated (p = 0.02) with increased CysLT receptor 1 (CysLTR1) methylation in peripheral blood mononuclear cells which, in turn, was marginally correlated (p = 0.06) with decreased CysLTR1 expression. Decreased CysLTR1 expression was, in turn, correlated (p = 0.0007) with FEV1 declines. During the same time period, increased methylation of GPR17 (a negative regulator of CysLTR1) was observed after DE exposure (p = 0.02); this methylation increase was correlated (p = 0.001) with decreased CysLTR1 methylation which, in turn, was marginally correlated (p = 0.06) with increased CysLTR1 expression; increased CysLTR1 expression was correlated (p = 0.0007) with FEV1 increases. Collectively, these data delineate a potential mechanistic pathway linking increased DE exposure-attributable CysLT levels to lung function declines through changes in CysLTR1-related methylation and gene expression.

The PedBE clock accurately estimates DNA methylation age in pediatric buccal cells.

The development of biological markers of aging has primarily focused on adult samples. Epigenetic clocks are a promising tool for measuring biological age that show impressive accuracy across most tissues and age ranges. In adults, deviations from the DNA methylation (DNAm) age prediction are correlated with several age-related phenotypes, such as mortality and frailty. In children, however, fewer such associations have been made, possibly because DNAm changes are more dynamic in pediatric populations as compared to adults. To address this gap, we aimed to develop a highly accurate, noninvasive, biological measure of age specific to pediatric samples using buccal epithelial cell DNAm. We gathered 1,721 genome-wide DNAm profiles from 11 different cohorts of typically developing individuals aged 0 to 20 y old. Elastic net penalized regression was used to select 94 CpG sites from a training dataset (n = 1,032), with performance assessed in a separate test dataset (n = 689). DNAm at these 94 CpG sites was highly predictive of age in the test cohort (median absolute error = 0.35 y). The Pediatric-Buccal-Epigenetic (PedBE) clock was characterized in additional cohorts, showcasing the accuracy in longitudinal data, the performance in nonbuccal tissues and adult age ranges, and the association with obstetric outcomes. The PedBE tool for measuring biological age in children might help in understanding the environmental and contextual factors that shape the DNA methylome during child development, and how it, in turn, might relate to child health and disease.

Systematic evaluation and validation of reference and library selection methods for deconvolution of cord blood DNA methylation data.

BACKGROUND: Umbilical cord blood (UCB) is commonly used in epigenome-wide association studies of prenatal exposures. Accounting for cell type composition is critical in such studies as it reduces confounding due to the cell specificity of DNA methylation (DNAm). In the absence of cell sorting information, statistical methods can be applied to deconvolve heterogeneous cell mixtures. Among these methods, reference-based approaches leverage age-appropriate cell-specific DNAm profiles to estimate cellular composition. In UCB, four reference datasets comprising DNAm signatures profiled in purified cell populations have been published using the Illumina 450 K and EPIC arrays. These datasets are biologically and technically different, and currently, there is no consensus on how to best apply them. Here, we systematically evaluate and compare these datasets and provide recommendations for reference-based UCB deconvolution. RESULTS: We first evaluated the four reference datasets to ascertain both the purity of the samples and the potential cell cross-contamination. We filtered samples and combined datasets to obtain a joint UCB reference. We selected deconvolution libraries using two different approaches: automatic selection using the top differentially methylated probes from the function pickCompProbes in minfi and a standardized library selected using the IDOL (Identifying Optimal Libraries) iterative algorithm. We compared the performance of each reference separately and in combination, using the two approaches for reference library selection, and validated the results in an independent cohort (Generation R Study, n = 191) with matched Fluorescence-Activated Cell Sorting measured cell counts. Strict filtering and combination of the references significantly improved the accuracy and efficiency of cell type estimates. Ultimately, the IDOL library outperformed the library from the automatic selection method implemented in pickCompProbes. CONCLUSION: These results have important implications for epigenetic studies in UCB as implementing this method will optimally reduce confounding due to cellular heterogeneity. This work provides guidelines for future reference-based UCB deconvolution and establishes a framework for combining reference datasets in other tissues.

Integrated analysis of environmental and genetic influences on cord blood DNA methylation in new-borns.

Epigenetic processes, including DNA methylation (DNAm), are among the mechanisms allowing integration of genetic and environmental factors to shape cellular function. While many studies have investigated either environmental or genetic contributions to DNAm, few have assessed their integrated effects. Here we examine the relative contributions of prenatal environmental factors and genotype on DNA methylation in neonatal blood at variably methylated regions (VMRs) in 4 independent cohorts (overall n = 2365). We use Akaike's information criterion to test which factors best explain variability of methylation in the cohort-specific VMRs: several prenatal environmental factors (E), genotypes in cis (G), or their additive (G + E) or interaction (GxE) effects. Genetic and environmental factors in combination best explain DNAm at the majority of VMRs. The CpGs best explained by either G, G + E or GxE are functionally distinct. The enrichment of genetic variants from GxE models in GWAS for complex disorders supports their importance for disease risk.

Glutaminase Deficiency Caused by Short Tandem Repeat Expansion in GLS.

We report an inborn error of metabolism caused by an expansion of a GCA-repeat tract in the 5' untranslated region of the gene encoding glutaminase (GLS) that was identified through detailed clinical and biochemical phenotyping, combined with whole-genome sequencing. The expansion was observed in three unrelated patients who presented with an early-onset delay in overall development, progressive ataxia, and elevated levels of glutamine. In addition to ataxia, one patient also showed cerebellar atrophy. The expansion was associated with a relative deficiency of GLS messenger RNA transcribed from the expanded allele, which probably resulted from repeat-mediated chromatin changes upstream of the GLS repeat. Our discovery underscores the importance of careful examination of regions of the genome that are typically excluded from or poorly captured by exome sequencing.

Genome-wide analysis of DNA methylation in relation to socioeconomic status during development and early adulthood.

OBJECTIVES: Socioeconomic status (SES) is a powerful determinant of health, but the underlying biological mechanisms are poorly understood. This study investigates whether levels of DNA methylation at CpG sites across the genome are associated with SES in a cohort of young adults in the Philippines. METHODS: DNA methylation was assayed with the Illumina HumanMethylation450 Bead Chip, in leukocytes from 489 participants in the Cebu Longitudinal Health and Nutrition Survey (mean age = 20.9 years). SES was measured in infancy/childhood and adulthood, and was based on composite measures of income, assets, and education. Genome-wide analysis of variable probes identified CpG sites significantly associated with SES after adjustment for multiple comparisons. Functional enrichment analysis was used to identify biological pathways associated with these sites. RESULTS: A total of 2,546 CpG sites, across 1,537 annotated genes, were differentially methylated in association with SES. In comparison with high SES, low SES was associated with increased methylation at 1,777 sites, and decreased methylation at 769 sites. Functional enrichment analysis identified over-representation of biological pathways related to immune function, skeletal development, and development of the nervous system. CONCLUSIONS: Socioeconomic status predicts DNA methylation at a large number of CpG sites across the genome. The scope of these associations is commensurate with the wide range of biological systems and health outcomes that are shaped by SES, and these findings suggest that DNA methylation may play an important role.

Systematic evaluation of DNA methylation age estimation with common preprocessing methods and the Infinium MethylationEPIC BeadChip array.

BACKGROUND: The capacity of technologies measuring DNA methylation (DNAm) is rapidly evolving, as are the options for applicable bioinformatics methods. The most commonly used DNAm microarray, the Illumina Infinium HumanMethylation450 (450K array), has recently been replaced by the Illumina Infinium HumanMethylationEPIC (EPIC array), nearly doubling the number of targeted CpG sites. Given that a subset of 450K CpG sites is absent on the EPIC array and that several tools for both data normalization and analyses were developed on the 450K array, it is important to assess their utility when applied to EPIC array data. One of the most commonly used 450K tools is the pan-tissue epigenetic clock, a multivariate predictor of biological age based on DNAm at 353 CpG sites. Of these CpGs, 19 are missing from the EPIC array, thus raising the question of whether EPIC data can be used to accurately estimate DNAm age. We also investigated a 71-CpG epigenetic age predictor, referred to as the Hannum method, which lacks 6 probes on the EPIC array. To evaluate these epigenetic clocks in EPIC data properly, a prior assessment of the effects of data preprocessing methods on DNAm age is also required. METHODS: DNAm was quantified, on both the 450K and EPIC platforms, from human primary monocytes derived from 172 individuals. We calculated DNAm age from raw, and three different preprocessed data forms to assess the effects of different processing methods on the DNAm age estimate. Using an additional cohort, we also investigated DNAm age of peripheral blood mononuclear cells, bronchoalveolar lavage, and bronchial brushing samples using the EPIC array. RESULTS: Using monocyte-derived data from subjects on both the 450K and EPIC, we found that DNAm age was highly correlated across both raw and preprocessing methods (r > 0.91). Thus, the correlation between chronological age and the DNAm age estimate is largely unaffected by platform differences and normalization methods. However, we found that the choice of normalization method and measurement platform can lead to a systematic offset in the age estimate which in turn leads to an increase in the median error. Comparing the 450K and EPIC DNAm age estimates, we observed that the median absolute difference was 1.44-3.10 years across preprocessing methods. CONCLUSIONS: Here, we have provided evidence that the epigenetic clock is resistant to the lack of 19 CpG sites missing from the EPIC array as well as highlighted the importance of considering the technical variance of the epigenetic when interpreting group differences below the reported error. Furthermore, our study highlights the utility of epigenetic age acceleration measure, the residuals from a linear regression of DNAm age on chronological age, as the resulting values are robust with respect to normalization methods and measurement platforms.